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We present a family of dynamic rotating cylindrically symmetric Ricci-flat grav- 
itational fields whose geodesic motions have the structure of gravitomagnetic jets. 
These correspond to helical motions of free test particles up and down parallel to the 



CJ , axis of cylindrical symmetry and are reminiscent of the motion of test charges in a 

^H ' magnetic field. The speed of a test particle in a gravitomagnetic jet asymptotically 

approaches the speed of light. Moreover, numerical evidence suggests that jets are 
^ . attractors. The possible implications of our results for the role of gravitomagnetism 

o : 

fSj in the formation of astrophysical jets are briefly discussed. 



PACS numbers: 04.20. Cv 



I. INTRODUCTION 



H \ A constant uniform magnetic field configuration in an inertial frame of reference has 

cylindrical symmetry; therefore, the motion of a test charge in this field is such that the 
particle's momentum and angular momentum in the direction of the field are constants of 
the motion. The particle in general moves with constant speed on a helix whose axis is along 
the field direction; moreover, the radius and step of the helical path are constants as well. 
The sense of helical motion about the direction of the magnetic field is positive (negative) 
for a test particle of negative (positive) electric charge. If the initial velocity of the particle 
is normal to the direction of the magnetic field, then the particle simply moves along a circle 
in the plane perpendicular to the field. 

The purpose of this paper is to study the analogous situation for the motion of free test 



particles in a gravitomagnetic field. The nonlinearity of this field implies that only a rough 
similarity may be anticipated. We investigate geodesic motion in a rotating dynamic space- 
time region that is a Ricci-fiat solution of Einstein's equations with cylindrical symmetry p|. 
Though this time- dependent gravitational case is considerably more complicated than the 
magnetic case, we find qualitatively similar phenomena. In fact, the helical motions up and 
down parallel to the axis of symmetry are reminiscent of the double-jet structure of certain 
high-energy astrophysical sources. 

Quasars and active galactic nuclei generally exhibit distinct relativistic outfiows. These 
jets are conjectured to originate from massive rotating black holes surrounded by accretion 
disks; the outfiows are focused beams of relativistic particles that proceed up and down 
along the rotation axis of the black hole (see, for instance, |2|). Similar phenomena have 
been observed in other high-energy astrophysical sources such as the Galactic X-ray binary 

iet formation, the more 



systems. While our results suggest a mechanism for astrophysica 

complicated physical process involves general relativistic MHD |3|, 

together with previous efforts p-llOj — contributes to the purely gravitational aspects of this 



j. The present work- 



fundamental problem in astrophysics. 

For a subclass of the Ricci-flat solutions under consideration, we show that the geodesic 
equations have families of special exact solutions that we call gravitomagnetic jets. More 
precisely, a gravitomagnetic jet is a set of special geodesies. These generally exhibit hehcal 
motions about the axis of cylindrical symmetry and their union is a non-compact connected 
invariant manifold that attracts all nearby geodesies. While we highlight features of these 
jets in the source-free cylindrical spacetime region of interest and provide strong (numerical) 
evidence that these families are attractors, the interesting question of the nature of the 
external matter currents that could generate such a gravitational field remains beyond the 
scope of our present investigation. 

According to general relativity, a rotating mass generates a relativistic, and hence non- 
Newtonian, gravitomagnetic field that is due to mass current. The exterior gravitomagnetic 



field of the Earth has recently been directly measured via Gravity Probe B (GP-B) |ll |. 



On the theoretical side, solutions of Einstein's equations in the case of axial symmetry have 



received attention for a long time (see Ch. VIII of 12| and references cited therein). In 



particular, rota ting solutions with cylindrical symmetry have been investigated by a number 



of authors (see 



13 



- Il5| and references therein). 



Previous interesting work on exact cylindrically symmetric gravitational fields in con- 
nection with the origin and structure of astrophysical jets has mainly involved the study of 



geodesies in the interior of time- independent rigidly rotating dust cylinders 16] . The behav- 
ior of free test particles in this case is directly influenced by the gravitational attraction of 
the rotating dust particles; to avoid this circumstance, we concentrate here on certain source- 
free gravitational fields that happen to depend exponentially upon time. These Ricci-fiat 
fields are perhaps more representative of the strongly time-dependent near-zone exteriors 
of accreting and growing gravitationally collapsed configurations where astrophysical jets 
are expected to originate. The temporal variation of the gravitational fields considered in 
our work leads to a significant and surprising feature of the special exact solutions of the 
geodesic equations. The speeds of test particles in gravitomagnetic jets start from values 
that are always above a certain minimum speed, rapidly increase along their paths and 
asymptotically approach the speed of light. The minimum speed represents the speed of 
circular motion perpendicular to the axis of rotation. The range of this minimum speed 
turns out to be from zero up to about 0.63 c. 

Is the main general relativistic problem of high-energy astrophysical jets solved in our 
paper? Our gravitational model is certainly too elementary to be an adequate representation 
of the complex physical situation. Nevertheless, in our simple model free test particles appear 
to be exponentially accelerated to almost the speed of light resulting in ultrarelativistic jet 
streams parallel and antiparallel to the axis of rotation. 

The plan of the paper is as follows. The class of spacetime metrics that we study in 
this paper is described in section II. Each metric in this class involves a function X{r) 
that is a solution of an ordinary differential equation and must be so chosen as to render 
our system of cylindrical coordinates admissible in the spacetime region of interest. The 
procedure for the determination of an appropriate X{r) is described in section III. The 
rotational aspects of the resulting gravitational field are discussed in section IV. Section V 
treats the motion of free test particles and null rays in the spacetime region of interest. This 
section also contains a discussion of the special analytic solutions of the geodesic equations 
that form gravitomagnetic jets. Numerical evidence that gravitomagnetic jets are indeed 
attractors is presented in section VI. Section VII contains a discussion of our results. For 
clarity of presentation, some of the detailed calculations and mathematical arguments are 
relegated to the appendices. Specifically, Appendix El contains useful formulas related to 



the principal dynamic spacetime metric under consideration in this paper. The proofs of 
the main mathematical results regarding the solutions of a certain highly nonlinear ordinary 
differential equation are given in Appendix[B], where, due to the nature of the subject matter, 
some of the notation employed is independent of the rest of the paper. Appendix O treats 
the geodesies of the time-reversed metric. 

II. SPACETIME METRIC 

In a study devoted to gravitational radiation [l|, a solution of the source-free gravi- 
tational field equations was described and partially interpreted using rotating cylindrical 
gravitational waves. The present paper is about the physical interpretation of a variant of 
this solution in a different physical domain. 

Consider the spacetime metric given in Eq. (35) of [1| 

_ ds^ = e-^—i-XHP + \dr'') + ^{IXdi + d<^f + e'dP (1) 

in (t, r, $, Z) coordinates. Here we employ the notation and conventions of Ref. [l|, so that 
the speed of light in vacuum is unity (c = 1) and the metric signature is +2. Moreover, 
Xr = dX/dr, £ is a constant and (r, $, Z) are standard circular cylindrical coordinates. The 
function X is a solution of the differential equation (cf. Eq. (31) of [l|) 

2 -ja2 dX . , 

dr"^ dr 

which is a special case of the generalized Emden- Fowler equation of the type y" = aq{x)y^y'"^ 
with y' = dy/dx, m = 1, n = —2, a = —1 and q{x) = x~^ (see [l7|). We transform it to the 
Lotka-Volterra system (JB27I) in Appendix [Bl We do not know an explicit general solution 
of Eq. ([2]). For the only explicit solutions known to exist, we note that X = constant is 
unacceptable as the 4D metric ([1]) would then degenerate into a 3D spacetime, and the 
solutions 

^ = ±(|'-)-"' (3) 



lil in |li 



correspond to the special free rotating gravitational waves investigated in detail in [ij and 
the references cited therein. 

About fifteen years ago, searching for a general description of rotating cylindrical grav- 
itational waves, one of us (BM) found the Ricci-flat metric ([1]), where X is a solution of 



Eq. ([2j) once i?^^ = 0. Among the solutions of Eq. ([2]), there is a class for which X^ is a 
monotonically decreasing function of r as well as a class where X^ is monotonically increas- 



ing. The 



fl 



ormer class of solutions, which includes Eq. ([3]) as a special case, was investigated 



in Ref. [l| and related to rotating gravitational waves. The physical interpretation of the 
latter class of solutions is taken up in the present work. 

Let us note that if X is a solution of Eq. ([2]), then so is —X; therefore, it is possible 
to assume £ > in Eq. ([1]) with no loss of generality. We define a lengthscale A such that 
i := A^^/^. Only positive square roots are considered throughout. Let us now introduce the 
dimensionless quantities f, X and z such that 

f = Ar, X = \-^'^X, ~z = \-^Z. (4) 

Moreover, t in Eq. ([1]) is dimensionless as well; therefore, we assume that the physical time 
coordinate is given by A't, where A' is in general a different arbitrary constant lengthscale. 
Under the scale transformation (r, X) i— )• (f, X), Eq. ([2]) remains invariant. Furthermore, 
with s = X^^s the spacetime interval ([T]) essentially remains invariant as well; that is, 
dropping all the tildes and working only with dimensionless quantities, the spacetime metric 
takes the form 

- ds^ = e-'^{-X^df + ^dr^) + —{Xdt + d<P? + e'dz\ (5) 

X r-^ r 

where t = t and $ = 0. Starting from dimensionless [t, r, 0, z) coordinates, we can always 
return to regular coordinates by choosing arbitrary lengthscales A and A'; then, the physical 
coordinates are (A't, A~^r, </>, A^;) and the other physical quantities in the metric are A^^^X 
and As. 

For the physical interpretation proposed in this work, the spacetime metric is obtained 
from Eq. ([5]) by t i— )■ — t. Henceforth, we will deal with dimensionless quantities and the 
spacetime metric 

- ds^ = e'^i-X^dt^ + \dr^) + -{-Xdt + d^f + e-'dz\ (6) 

where X is a solution of Eq. ([2]) specified in the next section, and we will show that under 
certain conditions geodesies of metric ([6]) allow gravitomagnetic jets. Thus the main focus 
of the following sections and Appendices lAl and IB] is on metric (|6]); we return to the time- 
reversed case — namely, metric ([5]) — in Appendix O 
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As shown in [l|, Eq. ([6]) represents an algebraically general Ricci-flat solution of type I 
in the Petrov classification. It admits two commuting spacelike Killing vector fields dz and 
d^ associated with the cylindrical symmetry of the gravitational field. The corresponding 
two-parameter isometry group is not orthogonally transitive. The invariant magnitude of 
the hyper surf ace-orthogonal dz is given by exp(— |t), while for d^, which is not hypersurface- 
orthogonal, the invariant magnitude is r~^/^ exp(it). Thus t and r can be invariantly defined 
in this way. The t = constant hypersurfaces are always spacelike. We interpret t as the time 
coordinate in this paper; therefore, at a given time t, p = r~^/^ is, up to a constant factor, 
an appropriate radial coordinate. Using this quantity, metric (^ takes the form 

- ds^ = -le'p^^i-X^dt^ + 4dp^) + e'p\-Xdt + d^f + e-'dz\ (7) 

so that the circumference of a circle orthogonal to the axis of cylindrical symmetry is 27rp*, 
where p* = pexp(|t). The symmetry axis is thus defined by p = (or r = oo). Despite its 
shortcoming as a proper radial coordinate, we nevertheless use r extensively in this paper 
to simplify formulas. 

We define the proper radial distance R* = exp(t/2)i? via Eq. ([7]) such that for the proper 
radial distance from the axis to a point with p = po, R = Rq is given by 

where tq = 1/po- 

The coordinates x'^ = (t, r, 0, z) are assumed to be admissible in the spacetime domain 
under consideration here. Thus, the gravitational potentials must be sufficiently smooth 
functions of these coordinates; moreover, the Lichnerowicz admissibility conditions 18| re- 
quire that the principal minors of the metric tensor and its inverse be negative for our choice 
of metric signature. For a symmetric n x n matrix M, the principal minors are given by 



det 



Mil ■■■ M; 



Ifc 



(9) 



Mki ■■■ Mkk 
for /c = 1, . . . , n. These admissibility conditions are essentially equivalent to the inequalities 

XirXr - X) > 0, XXr > 0, (10) 



for all values of r in the domain of the solution X. This fact follows from a detailed exam- 
ination of the spacetime metric tensor g^^, with —ds"^ = g^ydx^dx^ defined in display (|6]), 
and its inverse given in Appendix E] we note, in particular, that for an admissible solution 

—t 1 —t 

-9'' = ^>0, {-gr/' = '-^>0. (11) 

There are four algebraically independent scalar polynomial curvature invariants in a Ricci- 



flat spacetime that can be represented as (see Ch. 9 of |13[ |) 

Xl = R.upaR^"'" - tR.upaR*^"'", (12) 

X2 = R.^p^R^'^'^^R^p ^"^ + iRp,,„R'"'''^R\p ^r (13) 

For metric ([6]), Xi and X2 are both real and are given by 

Xi = 777-^(1 - 3rX2 - 2r^XXr - r^X'Xr + r^X^X^), (14) 

I, = - /I (I - 2rX^ - 2r'XXr + 2r'X'Xr + r^X'X^). (15) 

These invariants vanish in the case of special solution ([3]), but are in general nonzero. For 
any finite value of the coordinate time t, Xi and X2 depend on the particular branch of the 
function X. This will be discussed in detail in section III. 

III. CHOICE OF X 

The mathematical investigation needed to identify admissible solutions X{r) of Eq. ([2]) 
that would allow the existence of jets is given in Appendix [HI here, we present the main 
conclusions of this analysis. 

The character of metric tensor ([6]) depends on the choice of a solution of the differential 
Eq. ([2]). We will consider admissible solutions in (t,r, 0, 2) coordinates that correspond to 
solutions X of the differential Eq. ([2]) such that X'^ is monotonically increasing with r and 

Q{r) := rXr - X (16) 

is such that XQ > 0. 

Throughout this paper, we take advantage of the scale transformation already mentioned 
in the previous section: The differential Eq. ([2]) remains invariant under the scaling (r, X) h-)- 




FIG. 1: The top panel shows computer-generated plots of admissible solutions ±X versus r for the 
differential Eq. ([2]). The thick curve is for the initial data X{1) = 1 and Xr{l) = 2; the thin curve 
is for the data X{1) = 1 and Xj.(l) = 3. The middle (bottom) panel depicts computer generated 
graphs of F, H, W and X for the "thick" ("thin") solution. The graphs of F, H and X vanish 
at the left end-point r^ of the domain of definition of X, which is r^ ~ 0.720 for the thick solution 
that allows jets and r^, ~ 0.787 for the thin solution that does not allow jets. 
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{r,X), where for a > 0, f = ar and X{r) = a~^^'^X(r). Using scale-invariant variables, it is 
then possible to reduce the second-order Eq. ([2]) to equations of first order — see Appendix [Bi 
Indeed, consideration of scale-invariant quantities leads to substantial simplifications in our 
work. The admissible class that we seek is a scale-invariant subclass of solutions of Eq. ([2]) 
that are all of the form depicted for two cases in the top panel of Figure 1. In fact, in 
Appendix [B] we prove that there is an open set of initial conditions that corresponds to 
admissible solutions each of which exists on an interval (r?,, oo), where ri, > 0, with limiting 
values X{ri,) = and Xr{rh) = oo such that XXr = 1/r^ at rf,, the function X'^ is increasing, 
X^ is decreasing and X"^ approaches infinity as r approaches infinity. In addition, we show 
that there is an open set of initial conditions for solutions defined on (rf,, oo) with the same 
end-point asymptotics such that the scale-invariant function F defined by 

F{r)=r^X{r)Xr{r)-l (17) 

has a zero rj in the interior of the interval of existence (r^, oo). The physical significance of 
this property is that it is the necessary and sufficient condition for the existence of jets; that 
is, if F has such a zero, then the (timelike and null) geodesic equations admit special helical 
solutions propagating upward and downward on an interior cylinder of radius rj. These 
special solutions form gravitomagnetic jets. We also show that there is an open set of initial 
data corresponding to solutions that satisfy both conditions. 

To clarify the connection of the present work with the solutions treated in Ref. [ll], let us 
consider, as in [l|, the scale-invariant quantities 

A 2 AXr ,„. 



and note that these quantities are related by the autonomous differential equation 

da 3a{a — 6"^) 



(19) 



d6 2S{a - 5) ' 

By introducing a temporal parameter 6, we will also consider the corresponding first-order 
system 

f = 25ia-S), 

da (20) 

- = 3aia-n 

whose phase portrait is depicted in Fig. |2l With the choice of r as a positive radial coordinate 
in metric (E]), (5 > and hence the first and fourth quadrants of Fig. [2] are physically relevant. 
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a 




FIG. 2: Computer generated phase portrait for system (j20p . The rest point at {S,a) = (0,0) is a 
degenerate singularity of system ()20p . while the other rest point at (1, 1) is a simple singularity. 
This latter isolated singularity is a spiral point that corresponds to the special solution ([S]). 



The solutions in the first quadrant have a > or XXj. < 0, so that X"^ monotonically 
decreases with r. These solutions, which were discussed in [l|, tend to the special rest point 
{6, a) = (1, 1) that corresponds to exact solution ([3]). The present paper is therefore devoted 
to the study of solutions in the fourth quadrant for which a < and hence XXr > 0. Since 
system fl20p is autonomous and the coordinate axes are invariant, it follows immediately 
that the fourth quadrant is invariant; that is, solutions that start in the fourth quadrant 
stay there [l9|. In particular, X{r)Xr{r) > as long as such a solution exists. If X{r) > 0, 
then Xr{r) > 0, X increases as r increases, and Xrr{r) is negative. Thus, if X{r) > 
for some r > 0, this function increases as r increases and is concave down. Using our 
symmetry, if X{r) < for some r > 0, then X decreases with r and is concave up (see 
Fig. [T]). Each member of this class of solutions of Eq. (^ is asymptotic to the degenerate 
rest point {6, a) = (0, 0) that corresponds to the axis of cylindrical symmetry, where r = oo 
and X"^ = oo. In fact, it can be shown that, with our choice of temporal variable 6, each 
solution is asymptotic to (0, 0) (corresponding to axis of cylindrical symmetry) in positive 
time and approaches a point at infinity {6, a) = (oo, — oo) in negative time (corresponding 
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to the cylindrical boundary r = r^ > and XXr = r^^'^) 




1.5 



2.0 



2.5 



FIG. 3: Plot of R versus p corresponding to the spacetime defined by X with X{1) = 1 and 
Xr{l) = 2, as in the middle panel of Fig. [TJ 



Consider a hypersurface such that at every instant of time t, the hypersurface is char- 
acterized by the cylindrical surface r = constant for r;, < r < oo. The gravitomagnetic 
jets lie on such a hypersurface with r = rj. The normal to such a hypersurface is given by 
n^ = (0, 1,0,0); hence, n^n'^ = g'"" = exp{—t)r^X/Xr > using Eq. flA2p of Appendix A. 
That is, the hypersurface is always timelike, but as r — )■ r^^, n^ra'^ — > 0, so that the boundary 
hypersurface is null. 

The admissibility conditions require that we explicitly leave the inner and outer bound- 
aries out of the domain of the solution X. That is, the axis of cylindrical symmetry (r = oo) 
and the boundary cylinder (r = r^) are excluded from the spacetime region of physical 
interest, since (?** vanishes at the axis and {—g)^^''^ vanishes at the boundary (see Eq. ( ITTl) ). 

It is interesting to describe further some of the significant properties of the admissible 
solutions. Let us first note that the invariance of Eq. ([2]) under X i-> —X is a scale- 
invariant property; henceforth, we will work exclusively with the X > branch. We show in 
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Appendix [B] that near the axis (p — )■ and r — t- oo), 

XM~ar + aS,-J^ + ^ + 0(l). (21) 

where a > and Sq < are constants. This imphes, among other things, that the function 
Q > 0, which transforms hke X under scahng and satisfies dQ/dr = —Xr/{rX^), decreases 
monotonically from infinity at the boundary cyhnder to Q{oo) = —aSo > at the axis 
of symmetry. The behavior of X{r) near the boundary cyhnder is more comphcated (see 
Appendix [B]); we show there that for sufficiently small r — r^ > 0, 

r'XX^ = 1 + (-V^Air - nfl' + ^ii±^(r - r,) + 0{{r - nf'), (22) 

\rf,/ orb 

where A is a constant (cf. Eq. (1B44P ). This is a particularly useful relation as it involves 
the behavior of F{r) as r — > r^. For the admissible solutions that allow jets, A < 0. 

When A > 0, F{r) is positive for r in some interval whose left endpoint is Vb, while for 
A < 0, it is negative. The admissible solutions that allow jets are a subset of solutions of 
Eq. ([2]) for which A < 0. Using only the leading term of the series in Eq. fl22|) . we find the 
asymptotic differential equation 

Xir)Xr{r) = 1 (23) 

Its solution is given by 

X{r) = ±[X{rof-- + -Y^'. (24) 

r tq 

In the special case of the boundary cylinder where X{rb) = and r > r^, 

A-(r) = ±(!)-/^(^-l)"l (25) 

Based on this result, we let X{r) = e^/^x(e), where e = — 1 + r/vb; then, Eq. ([2]) implies 
a nonlinear second-order equation for x{^)- It turns out that this equation has a unique 
analytic solution near e = for A = 0, see Eq. ( JA12I) . For A ^ 0, there is no analytic 
solution; however, Eq. (1B43I) implies that as r — )■ r^^ 

X(,.)^±(i)"|e- + ^.+ <^.- + 0(.^),. (26) 

It follows from Eqs. (l2T])- (l22l) that for an admissible solution the proper radial distance 
from the axis to the boundary — namely, exp(t/2)i?(, — is finite at every given finite instant 
of time t. Moreover, near the axis p = 0, i? ~ p^, while R has infinite slope as p — )■ pb- 
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Figure [3] illustrates R versus p for the solution given in the middle panel of Fig. [1] We 
note that the requirement of elementary flatness is violated near the axis of cylindrical 
symmetry. That is, for an infinitesimal spacelike circle around the axis, the ratio of the 
circumference, 2-7rpexp(t/2), to the proper radius, p^exp(t/2), is 27r/p instead of 27r at a 
given instant of time t and diverges as p — )■ 0. For a thorough discussion of subtle issues 



regarding cylindrical symmetry and its axis, see J20|] and the references cited therein. We 
recall that in our analysis the axis (p = 0) is already excluded from the physical domain in 
order to satisfy the admissibility conditions. 

To help distinguish the admissible solutions that do allow jets from those that do not, 
it is useful to consider two other functions associated with X: H{r) and W{r). The scale- 
invariant H is defined by 

H{r) = rXQ -1 = F- rX^, (27) 

so that H{rh) = 0, while for r — )■ oo, H ~ —a^Sor"^; moreover, if X allows jets, — 1 < 
H{rj) < by Eq. (l27j) . The function W scales as X and is defined by 

/X \V2 
^ = Q-{^) ■ (28) 

Thus W is positive near the axis; in fact, W{oo) = —aSo > 0. For the value of W near the 
boundary, let us first note that W can be written as 

HQ-X 

and for r ^ r^ , HQ ^ A/rl^^ from Eq. ([22]). Thus W{n) = A/{2rl^^). The solutions 
that allow jets are among those for which A < 0. Figure [1] illustrates, via the behavior 
of these functions, the difference between the solutions that allow jets and those that do 
not. The top panel contains the graphs of two admissible solutions: The thick (thin) curve 
represents a solution that allows (does not allow) jets. The middle (bottom) panel illustrates 
the behavior of F, H and W for the admissible solution that allows (does not allow) jets. 

It is important to point out that for the admissible solutions, the curvature invariants 
Xi and X2 given by Eqs. ( IT4l) and ( TT5l) are indeed finite for finite values of time t and 
radial coordinate r G (r^, 00). We remark in passing that they diverge in the infinite past 
{t — ;> —00); however, this circumstance is consistent with the emergence of the universe from 
a singular state as in the standard models of cosmology. Using Eq. flT7|) . it is possible to 
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express these scalars as 

-2t 

Thus Xi(rj) = 4rj exp(— 2t) and X2(r j) = for any solution that allows jets. Specifically, for 
the case that allows jets depicted in the middle panel of Fig. [H Xi has one zero and X2 has 
two zeros in the interval (r^, 00), while for the case that does not allow jets depicted in the 
bottom panel of Fig. [1], Xi has one zero and X2 has no zeros in the interval (r?,, 00). Moreover, 
it is straightforward to show that as the boundary cylinder is approached (r — > r;,), 

Xi ^ (4 - A')e~''rl X, ^ -^e-^V,^ (31) 

while as the axis of cylindrical symmetry is approached (r -^ 00), 

Xi ^ e-2*5oa-2, X2 ^ -^e-3*a-^ (32) 

Once X has been properly chosen, we can turn to the treatment of gravitational physics 
in the corresponding singularity-free spacetime region. This is an open hollow cylindrical 
domain that expands; it has an inner boundary around the symmetry axis (r = 00) and an 
outer boundary (r = r(,). We therefore consider test particles and gyroscopes in the radial 
interval (r^, 00) in the rest of this paper. 

The rotational aspects of the spacetimes under consideration here provide the basis for 
the interesting features that free test particle motion can exhibit in these gravitational fields; 
therefore, we now turn to the gravitomagnetic properties of admissible solutions. 

IV. GRAVITOMAGNETIC FIELD 

ocal equivalence, in the linear approxima- 



The gravitational Larmor theorem implies a 
tion, between gravitomagnetism and rotation 21 



22( 1 . This correspondence may be employed 



in order to provide a definite measure of the gravitomagnetic field. It is therefore useful to 
study the precession of ideal test gyroscopes that are held at rest in the spacetime region of 
interest; the precession frequency may then be identified with the gravitomagnetic field for 
the class of observers that carry the gyros along their world lines. 
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To simplify matters, we consider the class of fundamental observers that are spatially at 
rest by definition; that is, x* is constant for each i = 1,2,3 for a fundamental observer in the 
physical spacetime region of interest. Therefore, the four- velocity field of the fundamental 
observers is given in (t, r, 0, z) coordinates by 

X\,^ = {{-gur'/',0,0,0). (33) 

These observers' natural orthonormal tetrad frame is given by A'^/^^ , where 

A'^(,) = (0,(^„)-l/^0,0), (34) 

X\^^ = is',0,stO), (35) 

A'^(,) = (0,0,0,e*/2) (36) 

are the corresponding spatial unit directions. Here, s* and s"^ are given by 

/ = -e-*/2(X,g)-i/2, s'f' = e-*/2(^)i/^ (37) 



^X, 



r 



The fundamental observers are accelerated; their acceleration tensor O0(a){i3) is defined via 

where T is the proper time along the world line x^{T) of a fundamental observer such that 
dx^/dT = A^/^N. In analogy with electrodynamics, the antisymmetric acceleration tensor 
consists of an "electric" part uj{t){i) = ^(i) and a "magnetic" part W(i)(j) = e(j)(j)(fc)r2'^'^\ The 
corresponding spacelike vectors are then A^ = A(^i)\n and fi'^ = ^2'^^'^A'^,^^. It follows from 
a detailed calculation that in the present case 

1 „w2/^^r\l/2X, -Xg2 



^'^^ 2 \ X J XXrQ ^ ' 

Ai^^^ = -\e-"\X^Q)-^l' (40) 

are the nonzero tetrad components of the translational acceleration A^ of the fundamental 
observers. These are finite in the physical region (r;,, oo); moreover, as the symmetry axis is 
approached 

Ar) ^ -^e-'/^S-o > 0, A<^) ^ -ie-*/2a(-5o)-^/' < 0, (41) 

while near the boundary r ^ r'l^ one can show, using the expressions given in the previous 
section, that 

^(r) ^ i^e-''^An, ^(0) ^ 0. (42) 
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(2XQ)- 




FIG. 4: Plots of {2XQ)~^ versus p = r~^'^ for the two cases of X given in the top panel of 
Fig. [TJ The upper curve here corresponds to initial data X{1) = 1 and ^^(l) = 2; the lower curve 
corresponds to X{1) = 1 and Xr{l) = 3. In general, the maximum of {2XQ)~^ occurs where 
H{r) = rXQ — 1 vanishes; this happens at r = 1 for the upper curve. 



Furthermore, Q'^ can be obtained from 



-t/2 



{xQr\ 



(43) 



which is the rotation frequency about the z axis of the local spatial frame of the fundamen- 
tal observers with respect to the local nonrotating — that is, Fermi- Walker transported — 
frame. All of the other nonzero components of the acceleration tensor can be obtained from 
Eqs. (l39])-g0]) and g3]). 

These considerations imply that ideal test gyroscopes carried along the world lines of the 
fundamental observers precess, with frequency —Q^^^ about the z axis, with respect to the 
natural spatial frame of the fundamental observers. Thus 

1 . 



fi^ = (0,0,0, 



2XQ' 



(44) 



characterizes the gravitomagnetic field in this case. It is variable in magnitude but constant 
in direction (parallel to the z axis). We note that {2XQ)^^ > vanishes along the axis 
of cylindrical symmetry and tends to rb/2 with infinite slope at the outer boundary of the 
region under consideration here (see Fig. Hj). The situation is more complicated, however, 
when we consider the gravitomagnetic components of the curvature tensor as measured by 
the fundamental observers. In particular, as demonstrated in Appendix \^ these do not all 
vanish along the symmetry axis. 

To illustrate these results further, let us consider a unit spacelike vector field V^ that is 
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carried by the fundamental observers and is orthogonal to the z axis; that is, 

V^ = co^^\\^^^m^\\^^. (45) 

If this is a gyro axis, then it is Fermi- Walker transported along x^(T), namely, 

-^ = {A.VnX\,y (46) 

It follows from a detailed calculation that Eq. (146!) is equivalent to the condition that 
dip/dT = —Q^''\ Thus the gyro rotates with proper frequency dip/dT about the z axis 
with respect to the spatial frame of the fundamental observers. 

The electric and magnetic components of the Riemann curvature tensor, projected on 
the tetrad frame of the fundamental observers, are given in Appendix [XI 

V. GEODESICS 

Spacetime geodesies are generally obtained from the condition that the spacetime interval 
along the geodesic path be an extremal, namely, 5/ ds = 0, where ds^ is given by Eq. (|6]). 
The geodesic equation then takes the form 

In this section we treat timelike and null geodesies in turn. 

A. Timelike Geodesies 

The geodesies of the dimensionless spacetime metric ([6]) depend on the choice of solution 
X of the differential Eq. ([2]) or the equivalent first-order system 

^ = F, ^— = -^. (48) 

dr ' dr r^X^ 

For an admissible X, we choose its positive branch X{r) > 0; moreover, to obtain the 
corresponding geodesic equations with respect to proper time r, we define U := dr/dr and 
note the reparameterization of system 



f-' f--' f-l^- («) 
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The components of the four-velocity vector of the free test particle along Killing vector 
fields are constants of geodesic motion. Therefore, due to cylindrical symmetry, there are 
two constants of the motion Cz = Qzadx^'/dr, the specific momentum in the z direction, and 
C^ = g^adx°'/dT, the specific angular momentum about the z axis; that is. 

By inserting these relations into the metric and simplifying, we find that 

or 
where 

y ■■= [^,u' + ^y + ^' + (^>^~")]""- (53) 

Our choice of positive sign in Eq. ( 152|) is in conformity with the notion that the temporal co- 
ordinate should monotonically increase with proper time along the world line of an observer. 
The Christoffel symbols Y'' that are needed for the radial geodesic equation — namely, the 
component of Eq. (1471) for the variation of the radial coordinate along the geodesic path — are 
given by 



^ tt - 


-lri-l^Q^§), 




-r -i 




p. rXQ 

4>t 2F ' 


V = 

rr 


1 + 3rX2 + r'^XY 


2r2X2 


T^r 


rX 


'-H - 


2Y' 



(54) 
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Thus a full set of differential equations for the timelike geodesies can be expressed as 

dr X ' 

^ = (7, 

dr 

dr 

"^^ C.e\ (55) 



dr 

dr 

dY UY 



-t\2 



dr r2X2 

In a constant magnetic field configuration, the motion of a test charge is a combination 
of uniform rectilinear motion along the field direction together with uniform circular motion 
around this direction. The analogous situation in the gravitational field under consideration 
is, however, much more complex. In particular, the free motion of a test particle purely 
parallel to the rotation axis is impossible in the physical region (r^, oo), since system (155|) 
does not have a solution for constant r and (j) coordinates. However, for a subset of admissible 
field configurations, special circular and helical motions are possible. 

In analogy with the magnetic case, let us look for geodesic motion that is confined to a 
cylinder of fixed radius r > 0; that is, we let f/ = and dU/dr = in system (l55l) . The 
latter equation, after division by {d^i/drY > 0, can be written as 

r«(|)^ + 2r;,(|)+r;,^o, (56) 

which is reminiscent of the quadratic equation usually encountered in discussions of the 



gravitomagnetic clock effect 23|. It is simple to show from Eqs. flM|l that 



and 



(r^^) r^^r^^ — — — (57) 



r:, = §[Q-{^r][QUjri (58) 



It follows from Eqs. (I56])-(I5HD that 

dt 



Q±(r/x)i/2- 



(59) 



20 

For an admissible solution X of Eq. (|2]), X{r) > 0, Y{r) > and Q{r) > 0; hence, the upper 
sign in Eq. f l59|) would always lead to helical motion in the negative sense about the z axis, 
while for the lower sign the sense of the motion depends on the sign of 

W{r):=Q-[-) . (60) 

This function has been discussed in Sec. Ill; it is given by A/{2r^ ) at r = r^ and asymptot- 
ically approaches —aSo > as r — )■ oo. In the class of admissible solutions, either W{r) > 
or W has a zero r^^, in the physical interval (r^, oo). This situation is illustrated in Fig. [TJ 
In the former case, the sense of helical motion would always be negative, while in the latter 
case we find from Eq. (1561) 

^rf0^-='-» " ~2Q' ^^^^ 

For r G {rh,rw), however, there could be helical motion in the positive sense. This is crucial 
since we find from the relations for dt/dr and dip/dr in system ( l55l) that 

dt _ V 

# ~X(y + Qre-*)' ^ ' 

In this relation we must have C,^ = 0, since it follows from Eq. fl59|l that dt/d(j) must be 
constant for fixed r. Thus, dt/d(j) = X~^ > and from Eq. (!59l) we find that 

~ Q - (y/x)V2 = X ^^^^ 

holds for rY = {Y/Xy/\ or 

r^XY = 1. (64) 

If this condition is satisfied for some rj G {rb,r^), then there is in general helical geodesic 
motion in the positive sense about the z axis regardless of the value of C^ 7^ 0. From the 
definition of F{r) in Eq. ( ITTll . we see that if F{r) = has a solution rj in the interior of the 
physical interval (r6,oo), then special helical solutions of the geodesic equation (H71) exist; 
these special solutions reside in the gravitomagnetic jet. 

For Cz = C^ = 0, there is a special family of timelike circular geodesic orbits at z = zq 
with radius rj given by Eq. (|6] 



= 0o + X(rj)(t-to) (65) 

and 

e*/' = e*°/2 + ^o(r-ro), (66) 
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where to, tq, 0o and zq are initial values of the corresponding quantities. These circular 
orbits separate the up and down helical motions about the z axis. Indeed, the special class 
of timelike geodesic orbits that define a gravitomagnetic jet include the special circular 
geodesies as limiting cases, for Cz -^ 0, of helical motions parallel {Cz > 0) and antiparallel 
{Cz < 0) to the axis of rotation. However, unless specified otherwise, such as in Sec. VI, 
for example, the circular geodesic orbits of constant speed constitute a relatively negligible 
set of measure zero and are therefore generally ignored in our discussion of gravitomagnetic 
jets. 

The existence of these special solutions of the timelike geodesic equation depends on 
whether there is a radial coordinate rj G (r;,, oo) for which r'jX{rj)Xr{rj) = 1. Given an 
admissible solution X — two examples of which are depicted in Fig. [1] — one can show that 
either there is only one such rj or there is none. To see this, we simply note that the function 
(r^Xr)"^ > starts from zero at r = r^ with infinite slope, has a maximum at 2rX^ = 1 
and then drops off to zero as r — > oo. Thus, X and (r^X^)"^ either do not intersect each 
other for r G [vb, oo) or they do so at exactly one point. 

By imposing the conditions that C^ = and r = rj, it follows immediately that r, U, 
X and Y remain at their initial values — in particular, the right-hand side of the differential 
equation for dU/dr vanishes — while 

^ = oX(o)(e-* + C,^)V2, (67) 

dr 
This system can be reduced to quadrature in elementary functions. For Cz 7^ 0, its flow is 
given by 

t = 21n(-^ sinh[^o|C,|(r - Tq) + sinh-^(|C,|e*°/2)]), 

= 0o + ^(o)(t-to), (68) 

z = zo + ^ r smh^il-r j\Cz\{r' - tq) + smh-\\Cz\e'^'^')] dr', 

where we consistently use zero subscripts to denote initial values. It is clear that t and 
increase as r — tq increases; however, z — zq increases for C^ > and decreases for Cz < 0, 
thereby leading to gravitomagnetic jets propagating up and down parallel to the rotation 
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axis. The function z{t) may be expressed as 

z{t) = zo + K+{t)-K+{to), (69) 



where 



j\Cz 



K+ = -±^{k+{1 + K+y - \n[K+ + (1 + K+Y']}, 



(70) 



B. Null Geodesies 

Let (^ be an affine parameter along the world line of a null geodesic. It follows from the 
existence of the spacelike Killing vectors d^ and d^p that 

where C*^ and C^ are constants of the motion. The spacetime path is null ( ds"^ = 0); hence, 

Let us define U and V^ such that U = dr/dC, and 

y--=[y,U' + ^}Cl + Clre~'^)\"\ (73) 

Then, the geodesic equations for a null path are 

d^_V_ 
dC~X' 



— — u^e , 



dU V^ UV V -^ 

rfr _ UY 



(74) 
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As before, with C^ = 0, it is possible to find an exact class of solutions of these equations 
once r'^XXr = 1 for some rj G (r^,, oo). Each null geodesic in this special class follows a 
helical trajectory in the positive sense on the cylindrical surface r = rj with Cz ^ and 

t-to = o|C'.|(C-Co), 

0-0o = oX(o)|C',|(C-Co), (75) 

For the example of X depicted in the middle panel of Fig. [H rj ^ 0.7739, X{rj) ^ 0.4288 
and Xr{rj) ^ 3.8941. Let us observe that z{t) for the special null geodesies coincides with 
the late-time behavior of special timelike geodesies with C^ 7^ 0, since K~^ is given, as t — )■ 00, 
by 

We remark that the special case of null circular geodesies is excluded here; that is, 
Cz = C^ = is not possible, since dt/dC, would then vanish and this is forbidden. 

A characteristic feature of the special (timelike and null) solutions of the geodesic equa- 
tions is that the sense of helical motion is always positive. This is due to our choice of the 
positive branch of solutions of Eq. ([2]), namely, X > 0. In fact, metric ([6]) remains invariant 
under X i— )■ —X and cf) i— )■ —cf). Thus if we work exclusively with the negative branch X < 
instead, the sense of helical motion will be negative. Hence the significant feature of gravit- 
omagnetic jets that must be emphasized here is simply that in the double-jet configuration, 
both jets have the same helical sense. Moreover, it follows from the solutions of the equa- 
tions of jet motion that there is a characteristic exponential dependence of \z — zo\ upon 
the azimuthal angle (j) — (J)q in gravitomagnetic jets. It is important to emphasize that the 
radius of the helical path of a gravitomagnetic jet is constant only in terms of r; in fact, the 
helix expands as its proper radius is given by exp(t/2)i?j. We note that helical motions in 



astrophysical jets have been the subject of recent investigations — see, for instance, 2J, |25| 
and references therein. 



Jets 



It is important to note that for t — )■ 00, the special helical timelike geodesies approach 
the special null geodesies asymptotically; in fact, this can be simply seen from the formal 
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correspondence between the respective geodesic equations. That is, systems (155|) and (^^ 
become formally equivalent once V and V in Eqs. f l53|) and f l73|) . respectively, take the same 
form; this actually happens when exp(— t) — )■ in Eq. ( 1531) . To gain physical perspective, 
however, let u^ = dx^ /dr be the four- velocity vector of a free test particle following a special 
timelike geodesic. With regard to the fundamental observers along the path of the particle, 

«'^ = «(")A^(^), (77) 

where u^"''^ = 7(1, v), v is the local velocity of the particle as measured by the fundamental 
observers and 7 is the corresponding Lorentz factor. Thus 7 = —u^X'^,^., which can be 
calculated using Eq. (!33|) and the result is 
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1 + Cfe' \ 1/2 



) , (78) 



.l + H{rj)J 
where H is given by Eq. fl271) . 

For a jet, Cz 7^ and hence the Lorentz factor for a jet is always larger than 7min! which 
is defined by Eq. ( !78l) with Cz = 0. Thus according to the fundamental observers, 7min is 
the Lorentz factor for a free test particle on a circular orbit of radius rj about the axis of 
cylindrical symmetry. Let us recall that for admissible solutions of Eq. (j2]), rXQ > and 
hence 1 + H{r) > for all r G (r^, 00). When jets are allowed in an admissible solution, 
there exists a unique rj G (r?,, 00) for which F{rj) = 0. From H = F — rX"^, we find that 
— 1 < H{rj) = —rjX'^{rj), so that 7niiii corresponds to a scale-invariant minimum speed 
/3min given by 

PL = rjX\rj). (79) 

Thus the jet speed is always greater than /^min- For the solution displayed in the middle 
panel of Figure [H /Smin ~ 0.3772. In general, the local velocity of the particle as determined 
by the fundamental observers can be written as v = (f^, v^, v^), where Vr = 0, v^ = /^min and 
■yVz = Cz exp(t/2). It is clear from Eq. (fTHl) that for C^ 7^ 0, 7 diverges as t — )■ 00, so that the 
local jet speed asymptotically approaches the speed of light according to the fundamental 
observers. This circumstance comes about due to the specific exponential dependence of 
metric ([H]) upon time t. This fascinating dynamical feature of the gravitational field is 
presumably caused by an exterior configuration whose characterization would necessitate a 
separate investigation. As a free test particle in a jet follows a helical path either up or down 
on a cylinder of radius rj, the gravitational potentials along its world line vary in time in just 
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such a way that the particle is apparently accelerated with its speed approaching the light 
speed for t — t- oo. For such extremely energetic test particles, however, our approximation 
scheme may break down at some point; that is, the gravitational influence of the test particle 
on the background spacetime may no longer be negligible. 

The preceding considerations may be employed to illustrate certain asymptotic charac- 
teristics of the gravitomagnetic jets using Eq. ( !75|) . Let z* be the proper distance parallel 
to the axis of cylindrical symmetry; then, dz* = exp(— 1/2) dz from Eq. ([6]). Moreover, an 
appropriate radial coordinate is p* = pexp(t/2). Therefore, along the special null geodesies 
we have 

^ = ^«"^ p--r-/'\"\ (80) 

It follows that 



\rf- (81) 

How can one numerically find admissible solutions of Eq. (|2]) that allow jets? For every 
solution with or without a jet, scaling provides a one-parameter family of solutions of the 
same kind. That is, every jet solution belongs to a one-parameter family of solutions due 
to the scaling property of Eq. (|2]); indeed, this differential equation remains invariant under 
the transformation r y-^ ar and X{r^ H- cr^^/^X(r) for o G (0, oo). The scaled solution 
is defined over the interval (err?,, oo); moreover, Q i— !■ Qja^l'^^ W H- W/a^^'^, H \-^ H and 
F h^ F. In particular, if a jet exists in the original solution at rj, the new scaled solution has 
a jet at arj. This property may thus be employed to set rj = 1 for every jet solution. One 
can then numerically integrate Eq. ([2]) with initial conditions such that rj = 1, X{1) = 1/d 
and Xr{l) = "(9, where "(9 > ^min- Here 'd~l^ is the maximum allowed value of /3min, which 
according to Fig. [7] of Appendix [B] is about 0.63; hence, ^min ~ 1-6. All such solutions are 
admissible according to the arguments presented in Appendix [Bl Moreover, Vb ^ 1 — ?9~^/2 
for 'i? ^ 1 in accordance with Eq. (125|) . As -i? — )■ oo, all (azimuthal) helical motions disappear 
and the special timelike and null geodesies become vertical. 

Finally, along the special timelike geodesic path with -u^ = dx'^/dr given by Eq. (1671) . 
consider an observer that carries an orthonormal parallel-propagated tetrad frame A^, ^ with 
^''(0) ~ ^^ ^^^ ^ spatial frame given by a set of unit gyro axes that can be expressed in 
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(t, r, 0, z) coordinates as 






AA* 
^^ (1) 


= (0, rfx(o)e-/^ cos 1 ryV*/^ sin 1 0), 


(82) 


^'m 


= (0, rf X(o)e-*/^in| -ry^e"*/^ cos | 0), 


(83) 


A^3) 


= (oC„ 0, oX(o)C„ e*(e-* + C2)i/2). 


(84) 



A free pointlike test gyroscope with spin S'^ carried by the observer along the path is then 
given by S^ = S^''^A^,^-,, where S^^\ i = 1,2,3, are constants. It is straightforward to verify 
that the requirements of orthogonality {Uf^S'^ = 0) and parallel transport {DS^ /dr = 0) 
are satisfied; moreover, S^S^ = S'(j)S''^*) is a constant of the motion. The spin vector in 
general undergoes damped precessional motion of frequency | with respect to time t; by 
contrast, the orbital frequency of geodesic motion is X{rj). The precessional motion decays 
exponentially; in fact, as t -> oo the special timelike geodesic with Cz ^ approaches a null 
geodesic and S^ ~ {S^^'^Cz/\Cz\)u^, as expected 26|. This follows from the fact that for 
t — !■ oo, A^,-^N and A^^g) asymptotically tend to zero and A^(-3^ ~ {Cz/\Cz\)u'^. The projection 
of the curvature tensor on A'^, ^ turns out to involve somewhat complicated functions of time 
t. In principle, the curvature coniponents as measured by the observer may be used to study 
the generalized Jacobi equation |6[ along the special geodesic world line, but that is beyond 
the scope of this paper. 

VI. NUMERICAL RESULTS: JETS ARE ATTRACTORS 

The special exact solutions of the geodesic equations that correspond to jets have vanish- 
ing (canonical) angular momentum C^ = and are confined to cylindrical surfaces of fixed 
r = rj] otherwise, they have arbitrary Cz, to, 4'o and zq for an admissible solution X that 
allows jets; for example, the one with X(l) = 1 and Xr{l) = 2 given in the middle panel 
of Fig. [H For this X, we have numerically integrated system (l55l) with Cz = ±1 and initial 
data to = 0, To = 1, f/o = 0, 00 = and 2;o = at To = 0; we find that from C^ = to 
C^ = 0.9, solutions are attracted to the jet while for C^ > 1, solutions are not attracted to 
the jet. The numerical results are presented in Figs. [5] and [61 

More generally, our numerical experiments suggest that the codimension-two submanifold 

{(t, r,0, 2, [/): r = rj and f/ = 0} (85) 
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FIG. 5: Schematic diagram illustrating the results of our numerical work. Right panel: For a range 
of parameters close to the jet parameters, the radial coordinate as a function of the particle's 
proper time starts from its initial value tq and reaches rj. It is then fixed at rj while the test 
particle executes helical motion up or down parallel to the rotation axis as in jets (see Fig. [6]). Left 
panel: Beyond a certain range in the canonical angular momentum C^, the test particle cannot be 
confined; that is, rj is bypassed and the particle soon reaches the boundary cylinder, thus leaving 
the spacetime region of interest. In constructing this figure, we have used the solution X given in 
the middle panel of Fig. [H therefore, rj, Ki 0.720, rj Ki 0.774 and ro = 1. 



is a gravitomagnetic jet. The flow on the invariant manifold (1851) has a simple geometric 
interpretation in physical space: each geodesic remains at a fixed radius rj from the axis 
of cylindrical symmetry. Solutions with Cz ^ spiral around the z axis on a helix that 
becomes unbounded as proper time approaches infinity. The measure-zero set of solutions 
with Cz = remains bounded in circular motion about the z axis. While we have not 
explored the entire parameter and state spaces, our numerical experiments confirm that 
the manifold f l85|) attracts all nearby geodesies. Thus, our experiments suggest that after 
transient motions nearby geodesies spiral about the z axis with radii approaching rj and, 
except for the negligible set with C^ = 0, become unbounded as proper time approaches 
infinity. 

More precisely, in our four- dimensional spacetime the geodesic flow takes place in an 
eight-dimensional state space (the tangent bundle of the spacetime). By the introduction of 
proper time, we consider only unit-speed geodesies. This reduces the state space to seven 
dimensions. Cylindrical symmetry implies that there are two integrals of the motion C^ 
and Cz] the state space is thereby reduced to five dimensions (t, r, 0, z, U). Since the jet is 
confined to the cylinder r = rj, it is a three-dimensional manifold parameterized by (t, 0, z). 
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FIG. 6: The result of integration of system (j55p for timelike geodesies attracted to jets. The initial 
data at To = are to = 0, ro = 1, Uq = 0, (/>o = 0, zq = 0, Xq = 1 and Yq = 2. The parameters 
are C,^ = 0.9 and Cz = ±1. The left-hand plot (jet going up) is for Cz = 1 and the right-hand 
plot (jet going down) is for Cz = —1. The coordinates are {x,y,z), where x = pcoscj), y = psiiKp 
and p = r~^'^; moreover, z = In | lnz| for the left-hand plot and z = — In | In |z|| for the right-hand 
plot. We use z instead of z for the sake of clarity. 



We also note that the third and fourth differential equations for and z in the geodesic 
equations ( !55|) may be decoupled from this system. Thus, the attraction properties of the 
manifold of special solutions correspond to the behavior of r and U. A geodesic is attracted 
to the manifold ( l85l) if r approaches rj and U approaches zero as r — )■ oo. 



29 



VII. DISCUSSION 

The cylindrical symmetry employed throughout this work has been a useful simplifying 
assumption. However, axial symmetry is expected to be a better approximation for the 
treatment of high-energy astrophysical jets that appear in circumstances involving a rotating 
collapsed configuration surrounded by a rotating accretion disk. On the other hand, in this 
case an analytical treatment appears to be prohibitively complicated. 

How could one generate the gravitational fields discussed in this paper? In electrody- 
namics, the magnetic field B inside an infinite circular cylinder is uniform and parallel to 
the axis of symmetry, provided the cylinder is surrounded by a uniform current sheet. The 
connection between the source and the interior field is given hj B = Airi/c, where i is the 
amount of electric current per unit length of the cylinder. This is a good approximation 
around the axis near the center of a long solenoid. In our gravitational case, we would expect 
that the source-free dynamic interior solution could be joined — perhaps along the inner and 
outer boundaries — to an exterior solution that could serve as the source of the interior field. 
However, finding such a source requires a separate investigation that is beyond the scope of 
this paper. In fact, it may be advantageous to look for gravitomagnetic jets in more realistic 
axisymmetric systems in general relativity. 

Despite these drawbacks, it is remarkable that general relativity permits the existence 
of Ricci-fiat rotating cylindrical gravitational fields that admit gravitomagnetic jets, which 
are formed from solutions of the geodesic equation corresponding to generally helical mo- 
tions of free test particles up and down parallel to the axis of symmetry with speeds that 
asymptotically approach the speed of light. Particle acceleration mechanisms are important 



in astrophysics 27|. We have shown that general relativity can in principle provide a purely 
gravitational mechanism for the directional acceleration of test particles to ultrarelativistic 
speeds. 
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Appendix A: Curvature Tensor 



For the Ricci-flat spacetime metric represented by Eq. (j6]) and x^ = (t, r, 0, z), the metric 
tensor and its inverse are given by 
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The nonzero components of the connection, modulo its symmetry (F" 
expressed as 
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(A4) 
(A5) 
(A6) 



together with the FJ^^, components given in Eq. ( l54l) of Sec. V. We note that Q and H have 
been defined in Eq. flT6|) and Eq. fl27j) . respectively. 

The components of the curvature tensor projected on the tetrad frame of fundamental 
observers, 

can be represented as a 6 x 6 matrix in the standard manner, 



n 



£ n 
n -s 



(A8) 



where S and "H are symmetric and traceless 3x3 matrices in a Ricci-flat spacetime. We iden- 
tify S and T-i respectively with the electric and magnetic components of spacetime curvature 
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according to the fundamental observers. The tidal matrix £ is given by 

Q + 3X -P 



£ 



AQXXr 



-P Q 

-2Q - 3X 



(A9) 



where P = {Q/{rX)y/^. 

The magnetic part of the curvature is given by 



AQ^XXr^ 
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(AlO) 



where V can be expressed, using H = rXQ — 1, as P = HP/X. 

Near the symmetry axis, r — )■ oo, we have X ~ ar, X,. ~ a and Q ~ —aSo, so that e^£ 
and e^H are in general nonzero constant matrices. Similarly, e^£ and e^Ti are in general 
nonzero constant matrices at the boundary cylinder as well. That is, near the boundary 
r — > Tfe, X -> 0, Xj. — )• cx) and F -^ 0. Moreover, H = F — rX^ and using Eq. ( 122|) . we find 
that 

Q = Ar'J'- (All) 

V (^ / r=rh 

We recall that A < when jets are present and A > in the absence of jets. In the special 



case of y4 = 0, e'^'H vanishes at the boundary; moreover, it is possible to show that for r 



n, 



^(f)"(^4-l'^-«(''')- 



X = ±i-r^(l + ^e- ^e' + 0(e')). (A12) 

where r = rfe(l + e). This result corrects an error in Eq. (34) of [l|, where 3/76 occurs in 
place of 5/72. 



Appendix B: Properties of solutions of r X Xrr + X,, 







To avoid confusion, we emphasize that some of the notation employed in this appendix 
is specific to the mathematical arguments at hand and is independent of the rest of the 
appendixes or this paper. 

We begin with two obvious facts about the solutions of the differential Eq. ([2]): If X is a 
solution, then —X is a solution. Also, we have that 



dr 



{xD 



-& 



(Bl) 
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Given tq > and initial data X{rQ) > and Xr^ro) > 0, there is, by the usual existence 
theory, a unique solution of the differential Eq. ([2]) defined on some maximal interval {vf,, tb) 
with Tb < To < vb- In fact, tb = oo. To prove this, we may use the extension theorem 



for ordinary differential equations (see, for example, [19|). In effect, a solution continues 
to exist until it reaches the domain of definition of the differential equation or it blows up 
to infinity. By Eq. ( jBip and for our choice of (positive) initial data, X^. is a monotonically 
decreasing function of r as long as X{r) and Xr{r) are both positive. If either of these 
functions vanishes for some r > ro, then there must be a point where X^ vanishes. Let r^ 
be the infimum of all such points and note that by continuity Xr{r^^) = 0. The function 
X{r) = X{r^) is clearly a solution of the differential equation that has the same (initial) 
data at r^. By uniqueness the two solutions must be the same. But, this is a contradiction: 
the original solution is not constant. This proves that Xj. is a monotonically decreasing 
function of r on the entire interval of existence. In particular, we have that 

Xr{r) < Xriro) 

on this interval; therefore, 

X(r)<X(ro)+X,(ro)(r-ro) 

for all r > Tq. That is, X does not blow up at some finite r. Since X and Xr do not blow 
up for finite r, the solution may be extended to the interval (r^,, oo). 

We claim that rf, > and as r — ;■ r^ the function X approaches zero and X^ approaches 
infinity. Under our assumption that X and Xr are positive at the initial point ro, X is 
decreasing and Xr is increasing as r decreases toward the left-hand endpoint r^ of its maximal 
interval of existence. Suppose there is some point p in the interior of this interval such that 
X{p) = 0. The function X is then defined in an open interval containing p; therefore, Xrr{p) 
is finite. By inspection of the differential equation we must have Xr{p) = 0. But, this is 
impossible because Xr increases from its positive initial value. Thus, we may assume that 
X > on its maximal interval of existence. If r?, < 0, then r = is an interior point of the 
interval of existence and again Xr{0) = 0, in contradiction. Therefore, rf, > and X > 
on its maximal interval of existence. Suppose that X is bounded above zero on this interval 
and rj, > 0. Then, in view of equation ( jBip . for r^ < r < Tq and 

b := inf -r^X^{r), 

rt<r<ro 2 
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we have the inequahty 

-i:(^rir))<lx^{r) (B2) 

or 

|-ln(X,^(r))>-l. (B3) 

By integrating both sides of this inequahty on the interval (r, tq) and rearranging, we find 
that 

X2(r) < e-(^-^^«)/^X2(ro). (B4) 

In particular, Xr is bounded on the interval (r6,ro). But, since the interval of existence 
was chosen to be maximal, this is impossible by the standard extension theorem for ordi- 
nary differential equations: solutions continue to exist until they reach the boundary of the 
domain of definition of the differential equation or they become unbounded. In effect, the 
solution under consideration here must continue to exist unless X reaches zero, which we 
have excluded, or X^ grows to infinity. But, using our assumption that X is bounded above 
zero, we have proved that Xr is bounded; thus, we have reached a contradiction. It follows 
that X approaches zero as r approaches rf,. 

We will show that X^. blows up to infinity as r approaches r^. By writing the differential 
equation in the form 

integrating both sides from r to ro and rearranging, we find that 

Xr{r) = Xr{ro) exp ( ^ ° ^^^ ds) . (B6) 

Because Xrr < 0, the graph of X is concave down; therefore, this graph lies below each of its 
tangent lines. We know that Xr is increasing as r — )■ r^. If Xr is unbounded on the interval 
(rfe,ro), we have \im.^^^.+ Xr{r) = oo, as desired. On the other hand, if Xr is bounded then 
lim.^^^+ Xr{r) = K := sup^^^^^^^^ ^^(r) < oo. It follows that X{r) < K{r — Vh), where 
r i-> K{r — Tfe) is the function giving the limit tangent line to X at r;,. By inserting the 
inequality 



X2(r) K^{r-ny 
in Eq. (1B6I) . it is easy to see that 



AV(r)>X.Mexp(^y^ j— -jj d,) (B8) 



34 

for Th < r < vq. But, the last inequality implies that Xr{r) — )■ oo as r — )> r^, in contradiction 
to the assumption that X^ is bounded on (r(,,ro). 

To determine the admissible solutions that allow jets, we consider the new scale-invariant 
functions 

X = r^^^Xrir), y = r^^'^X{r) (B9) 



(first found by Weishi Liu 28|]) and note that 



dx X ,3 1 . dy 1 , 1 . ,^ . 

:r = -9--' / = -0^ + o^- BIO 

dr r 2 y^ dr r 2 

Thus, we have determined a first-order system equivalent to the second-order differential 
Eq. ([2]). By the change of independent variable r = e^, this system is made autonomous. In 
fact, with 

as) = x{e'), r^{s)=y{e'), (Bll) 

system (IBlOp is transformed to the autonomous system 



?7^ 2 2 



e = :T(o^'-i)' V = ^ + 7;V, (B12) 



where the overdot denotes differentiation with respect to the new independent variable s. 

Recall that solutions are admissible if XXr > and X{rXr — X) > 0. We will consider 
only solutions with X > for simplicity. Solutions with X < have similar properties by 
symmetry. The first condition is satisfied if both ^ and rj are positive; that is, solutions that 
start in the (open) first quadrant of the [C,, t]) space remain there. Because, 

X{rXr-X) = y-{x-y), (B13) 

r 

a positive solution of the differential Eq. ([2]) is admissible exactly when it is confined to the 
sector S in the open first quadrant where C, > f]- 

The positive ^ axis is part of the boundary of the domain of definition of system ( JB12I) . 
In particular, solutions of this system starting in S do not exit along this ray. An easy 
calculation using system f lB12p yields the inequality 

(e-^)|^=, = -^- (B14) 

It follows that S is not positively invariant because solutions that meet the upper boundary 
of S leave this sector. 
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An admissible solution (that is, one that stays in S) allows jets exactly when F{r) = 
r'^X{r)Xr{r) — 1 vanishes at some r > r^. In the new coordinates, 

r^XXr-l = xy-l. (B15) 



From Eq. flB12p . we have that 



-f(er/-l) =2>0. (B16) 

as €'7=1 

Therefore, solutions of system (1B12I) starting in the sector S cross the hyperbola .^r/ = 1 at 

most once. 

It is convenient to consider the change of coordinates 

u=^, ^ = |, (BIT) 



where the line u = may be viewed as the line at infinity (see 29|). Using the new 
coordinates, system (]B12p is transformed to 



u = -^{-v^ - n^), {; = -{v -v^ + u^). (B18) 

It has the same phase portrait in the first quadrant as the system 

3 

u = u{u^ - -v"^), v = v{y -v"^ + u^). (B19) 

The upper boundary of S corresponds to the line w = 1 in the new coordinates, and the 
curve xy = 1 (representing F{r) = 0) corresponds to f = u"^. 

We note that system ( JB19I) has two rest points on the line at infinity: (m, v) = (0, 0) 
corresponding to the boundary cylinder at r = r^, and {u,v) = (0, 1) corresponding to the 
axis of symmetry at r = oo. The second rest point is a hyperbolic sink, and the system 
matrix of its linearization is diagonal. The eigenvalue of the system matrix corresponding to 
the vertical direction is —3/2 and the eigenvalue corresponding to the horizontal direction 
is —1. It follows from basic invariant manifold theory that there is an analytic solution Z 
which approaches the rest point (0, 1) tangent to the line v = 1 and has the Taylor series 
expansion 

Therefore, Z approaches the rest point (0, 1) from below the line v = 1. Numerical experi- 
ments suggest that Z also approaches the origin (in the backward direction of the indepen- 
dent variable) and it crosses the curve v = v? (see Fig. [7]). 
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V 0.6 




FIG. 7: A portion of the {u,v) state space for system (JBIOP is shown here. We have plotted the 
hne V = 1, the parabola v = u^, the line v = (2/3) ^'^n (which is the vertical isocline) and the thick 
curve Z that is an approximation of the solution that approaches the rest point at (0,1) tangent 
to f = 1. This solution Z connects the two rest points and crosses the parabola. 



Our numerical evidence suggests that all solutions of system ( 1B19I) starting in the first 
quadrant to the left of the curve Z correspond to admissible solutions that allow jets. It is 
clear that every solution starting to the left of Z remains to the left; therefore, such a solution 
does not cross the line v = 1. This means that corresponding solutions remain in S for all 
time. As previously mentioned, this fact implies that the quantity X{r){rXr{r) — X{r)) 
is positive for all r so that these solutions are admissible. Figure [7] strongly suggests that 
an open segment of the parabola v = u^ lies to the left of Z. It follows that all solutions 
starting on this portion of the parabola correspond to spacetimes that allow jets. 

We note that Fig. [7] indicates the range of /Smin defined in Eq. fl79|) . Indeed, using the 
coordinates defined in this appendix. 



/3m 



y = r] 



V 

u 



(B20) 
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at points in {u,v) space on the curve v = u"^, which corresponds to evaluation at r = rj. 
Thus, we have that 

/^Min = u (B21) 

for those values of u such that {u, u^) is in the region of admissible solutions to the left of 
the curve Z depicted in Fig. [71 The corresponding range of u is approximately the interval 
(0,0.63). 

The nonlinear system (IB19|) will behave asymptotically the same as its linearization near 
the hyperbolic rest point (0, 1), which is given in the fourth quadrant by 

u = -\u, v = -v. 

Its solutions lie on curves of the form 

V = -a^f/^/^; 

hence, the corresponding solutions of the nonlinear system lie on curves with asymptotic 
expansions of the form 

Also, since the solution Z is analytic and tangent to the horizontal axis, a computation with 
power series can be used to show that it lies on an analytic curve of the form 

v = l- -u^ + 0{u^) 

near this rest point. Thus, all the admissible solutions have expansions of the form 

For t; = 1, which corresponds to the leading-order behavior of an admissible solution as 
r — )■ oo, we have that y = x and X{r) = rXr{r). The solutions of this differential equation 
have the form X(r) = ar, for some constant a > 0. This suggests the asymptotic behavior 
of the differential Eq. ([2]) as r — )■ oo might be obtained by studying solutions of the form 
X{r) = arS{\/r) for some function S such that S'(O) = 1. By inserting this relation into 
the original differential Eq. ([2]) and using rj = 1/r, we find that if 5*0 is a constant and S 
satisfies the initial value problem 

a^S^S" -7]^S' + 7]S = 0, 5(0) = 1, S'{0) = So, (B22) 
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then X{r) = arS{\/r) is a solution of (|2]). We note that the differential equation is not 
singular at S* = 1. Therefore the initial value problem has an analytic solution, which may 
be obtained by inserting a formal series representation for S in powers of tj and equating 
coefficients. The solution has the form 



^(r/) = 1 - -^ + 0{^') + Soilil + ^ + 0{r^')). 



It corresponds to 



X[r) = ar{l + ^ - + ^ + 0(-)). (B23) 



We have that limr_^oo''"-^r(^) — ^('") = —aSo. But for admissible solutions, the quantity 
rXr{r) — X{r) has the same sign as X{r); thus, we see that —aSo must be chosen to have 
the sign of X{r). That is, under our assumption that X{r) is positive, we take 5*0 < 0. 

While the numerical evidence is compelling, we will prove that there is an open set of 
admissible spacetimes that allow jets. Our argument uses new coordinates. We note first 
that for the basic differential Eq. ([2]), we have X(rf,) = and Xr{rb) = oo. To avoid 
the infinite derivative, we consider instead the inverse function Y, which is defined by the 
relations Y{X{r)) = r and X{Y{s)) = s, and note that it satisfies the initial value problem 

s^Y^Yss = Y,\ Y{0) = n, Ys{0) = 0. (B24) 

We define the scale-invariant quantities 

x{s) = s^Ysis), y{s) = s^Y{s) 

and 

e(t) = x(e-*), r/(t)=y(e-*) 

to obtain the corresponding first-order system 



e 



e = -^(3r/' + 0, r) = -(e + 2r/). (B25) 



It is convenient to make the change of coordinates 



q q 

that transforms system flB25l) to 



e = -, V = - (B26) 



p = p{p- q - I), q = q{l + 2p- q). (B27) 
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We note that system (1B27P is a Lotka-Volterra system, a family of differential equations 



that arises in many applications such as, for instance, chemical kinetics and ecology 30 1. 

Recall that the admissibility conditions are XXr > and X{rXr — X) > 0. The first 
condition is effectively redundant. Indeed the second condition is equivalent to 

r^^-l>0 
X(r) 

which can only be satisfied for XXj. > 0. Note that 

Xr{r) Y{s) y{s) r]{t) 1 



X{r) sY,{s) x{s) at) P{ty 
hence 



pW = ^.' «W = r^- 'S2«> 



Lf. 



where t = — InX. Thus, a solution is admissible when < p{t) < 1 for all t. Similarly, 
an admissible solution allows jets when r'^X{r)Xr{r) — 1 vanishes for some r > r^. This 
condition is satisfied for those solutions such that q{t) = 1 for some finite t, where t = — cxo 
at the symmetry axis and t = oo at the boundary. 

To show that there are admissible solutions that allow jets, we analyze the phase portrait 
of system (IB27P — see Fig.|Hl There are several important properties: (i) There is a hyperbolic 
saddle point at the origin whose stable manifold is the p axis and whose unstable manifold 
is the q axis. In particular, the axes are invariant sets, {ii) There is a hyperbolic source at 
{p,q) = (1,0), which corresponds to the axis of symmetry (r = oo), and a smooth solution 
that approaches this rest point in negative time tangent to the line q = — 2p + 2 with 
both p and q components positive, (iii) There is a hyperbolic sink at {p,q) = (0, 1), which 
corresponds to the boundary (r = rj,), and a smooth solution that approaches this rest point 
in positive time tangent to the line q = —2p + 1 with both p and q components positive. 
Moreover, this solution approaches the sink faster than the solutions along the q axis, (iv) 
Solutions cross the segment of the line p + q = 1 between this source and sink from below to 
above. More precisely, they cross in the positive direction of the normal to the line, which 
has components (1,1). {v) Solutions cross the line p = 1 from right to left, (vi) Solutions 
cross the line q = 1 from below to above. Facts (i), (ii) and {Hi) are proved by linearization 
and stable manifold theory. Facts (iv), (v) and (vi) are proved by checking the direction of 
the vector field associated with system f lB27|) on the lines. 
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FIG. 8: A portion of the phase portrait of system ()B27p is shown. The shaded region is bounded 
by the trajectory J that leaves the source at (p, q) = (1,0) tangent to the hne q = — 2p+2 (which is 
also depicted) and the trajectory that approaches the sink at (0, 1) tangent to the line q = —2p + 1 
(which is also depicted). Solutions that correspond to trajectories that lie in the shaded region are 
admissible (because their p components remain in the interval < p < 1) and allow jets (because 
they cross the line q = 1). 



Given the properties just enumerated, consider a solution J' leaving the source tangent 
to the line q = —2p + 2 and note that JT" starts above the line p + q = 1. Hence, by (iv) 
it stays above the line. Also by (v) and (vi), J is confined to the triangle with corners 
(1,0), (1, 1) and (0, 1) unless it crosses the line q = \. It is easy to see — by inspection or 
one may use the Poincare-Bendixson theorem — that J cannot remain in the triangle for all 
time unless it approaches the sink. In fact, it does approach this rest point but not until 
after it crosses the line q = \. To prove this fact requires a further analysis of the local 
behavior at the sink. The important point is that J approaches the sink above the line 
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q = —2p +1. A solution of the system approaching this hyperbolic rest point tangent to 
this hne — according to {in) — approaches faster than the solutions approaching along the 
q axis. It follows that all solutions except those on the trajectory approaching tangent to 
q = —2p + l must approach the sink tangent to the q axis. Hence, the solutions approaching 
above this tangent must cross the line g = 1, as required. 

It is interesting to note here the connection between these results and Eq. (!22|) . which 
implicitly expresses the way X{r) approaches zero as r — )■ r^ relative to the parameter A. The 
trajectory that approaches the sink at the boundary (p, q) = (0, 1) tangent to g = —2p + 1 
uniquely corresponds to A = 0, see Eq. ( JA12I) . All other trajectories must approach this sink 
tangent to the q axis and have A ^ 0; while those with A > stay below q = 1 and have no 
jets, the admissible solutions that allow jets all have A < 0. As shown in Fig. [HI there is a 
negative critical value of A (corresponding to trajectory J') below which the corresponding 
solutions are not admissible. 

If JT" crosses the line g = 1 at a point with coordinates (j, 1), then all solutions crossing 
the line g = 1 with p coordinate in the interval (0, j) are admissible and allow jets. It follows 
that j is the same as the maximum possible /S^in, which is approximately 0.4 according to 
Figs. El and El 

Let us return to the behavior of the solutions of Eq. ([2]) as r — )■ r^ to prove Eq. (1221) . 
This series representation may be found formally in several ways. To prove that it is correct 
seems to require several steps that are outlined here. 

In the coordinates of system ( JB27I) . which are related to solutions of Eq. ([2]) in dis- 
play (1B28I) . the corresponding problem is the asymptotic behavior near the rest point at 
{p, q) = (0, 1). For computational convenience — that is, the ability to compute at the origin 
with a system whose linear part is diagonal — we change coordinates via 

P = p, Q = q-l + 2p 

to obtain the equivalent system 

P = -2P + 3P2 -PQ, Q = -Q- 2P2 + APQ - Q^. (B29) 

Since the eigenvalue —2 of the coefficient matrix of the linear part of this system is twice 



the other eigenvalue —1, normal form theory (see 3l|) implies that there is a near-identity 



analytic transformation defined in a neighborhood of the origin that transforms system flB29p 
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to 

Xi = — 2xi + 7X2, X2 = —X2 (B30) 

for some constant 7. By performing the reduction to normal form to second order, we find 
that 7 = for system (1B29I) . More precisely, the near-identity change of coordinates 

xi = P - ^P' + PQ, x2 = Q + Ip^- 2PQ + g2 

transforms system (1B29P to a system with the same linear part and a nonlinear part without 
quadratic terms. Stated differently, our system is linearizable in a neighborhood of the origin; 
that is, there exists a near-identity analytic change of coordinates defined in a neighborhood 
of the origin that transforms system (1B29|) to (]B30P with 7 = 0. This fact implies that 
system flB29l) has an analytic first integral defined in a neighborhood of the origin. More 
precisely, if xi = P+G{P, Q) and X2 = Q+H{P, Q) is the near-identity change of coordinates 
(where G and H contain no linear terms), we may simply transform the analytic integral 
X2 — kxi = for system ( 1B30I) with 7 = to 

iQ + HiP,Q)f = kiP + GiP,Q)) 

to obtain the desired analytic first integral of system ( 1B29I) . Since the physical domain 
corresponds to P > 0, the corresponding constant k is positive. We may therefore define 
k = A^, where A can be positive or negative. Using this fact together with the implicit 
function theorem, it follows that the differential equation 

dQ _ -Q- 2P^ + 4PQ - g^ 
dP ~ -2P + 3P2 - PQ 

has an analytic family of solutions passing through the singularity at the origin. In fact, 



these solutions have the form 



^ 'n2 ^3 , 2^ +5^4 , r^fr^5\ 



P = -^{Q' -Q' + -i^^Q + om)- (B31) 

We will determine the behavior of X as r — > r^,. To this end, we assume that A ^ 0, 
view the integral ( 1B31I) in the abstract form P = f{Q, A^) and change coordinates back to 
{p, g)-coordinates with ( := q — 1 to obtain the relation 

p-fiC + 2p,A')=0. (B32) 
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Let J^(p,C, A) = P- fi( + 2p, ^2). We have J^(0,0, A) = and Jp(0,0, A) = 1; hence, 
by the iniphcit function theorem, we may solve for p as an analytic function of ( and A at 
(p, C, A) = (0, 0, A). Let us write the result as 

p = g{C,A'). (B33) 

This relation must hold for solutions of system (1B27P near the rest point at {p,q) = (0, 1). 
Using the second differential equation in system (1B27I) . we have that 



C = (C + 1)(2(7(C, A') -() = -( + OiC). (B34) 

The right-hand side of this differential equation is an analytic function of ( defined in a 
neighborhood of C = 0. Of course, we may recover solutions of system ( 1B27I) simply as 

p{t)=giat),A'), g(t) = l + C(t). 

Using the coordinates defined in displays (IB24I) and (1B28I) and noting that s = e~*, we 
have 

(B35) 



X(r) = e-, rX^(r) = ^ = i^^' 



q{t) i + CW ■ 

The linearization of the scalar differential Eq. (IB34P has eigenvalue —1; therefore, this 
differential equation is linearizable. More precisely, there is a near-identity analytic change 
of coordinates 

C = (t) + H{ct)) (B36) 

(defined near C = 0) that transforms Eq. flB34p to = —0, which has the general solution 



0(t) = 0o6xp(— t), where 0o = 0(0) is an integration constant. Thus, Eq. (1B34P has the 
corresponding solution 

C(t) = 0oe-* + i/(</)oe-*). 

In view of Eq. f lB35l) , we have 

We will show that it is possible to solve for X{r) in this equation. To this end, we first 
determine the series representation of H at the origin. Using Eq. (JB36I) . we have the relation 



C = (1 + H^{(P, A'))<f^ = -0(1 + H^{<p, A')). 
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After substitution for ( (again using Eq. (1B34P ). we may determine the coefficients of the 
power series representation of H. Using the power series in Eq. (IB37p . we find the relation 



which may be simphfied to 

Since X(rb) = 0, it follows that 0q = A^ri,. Thus, 0o = ^Avf^ ; however, for the sake of 
consistency with Eq. (1221) . we choose (pQ = —Ar^ . Thus, Eq. flB38p becomes 

x^ (y - ^^) + ^(^') = ^-r^- (B39) 



To solve for X(r) in Eq. flB39p . define 



r2 Ar^l'^ 



^(X,r)=r-r,-X2(^-^x)+0(X 



and note that 

Under the assumption that r^ 7^ 0, a condition that holds for admissible solutions that allow 
jets, the Weierstrass preparation theorem (see, for example, [19|) implies that 

Q{X, r) = {a{r) + h{r)X + X'')U{X, r), (B40) 

where the functions a, h and U are analytic in a neighborhood of X = and r = r^ and 
U{0,rh) 7^ 0. It follows that the left-hand side of Eq. (1B40I) is zero for {X,r) near (0,rf,) 
exactly when a + bX + X^ = or 

X = -(-6(r) ± {b^{r) - 4a(r))^/2). (B41) 

We return to Eq. f lB40p . expand both sides in power series of both variables about the 
point X = and r = ri, and equate coefficients to get 



2, , 2(9-^2^ 
— r - n) -^ 



(^i^) = — 2 (^ ~ ^b) ^71 {r - n) + 0{{r - n) ) 

4:A 4:A{9-A 



(B42) 
Kr) = ^(r - n) - ^^\^" ,;; \ r - nY + 0{{r - rfe)=^). 



3r^ 



3/2 ^ "^ 135r5/2 
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Substitution of these results in Eq. (lB4ip results in 



21/2 2 A . . 1(9 + A 

2 A(9-A2 



^W = ^(r - r.)V^ + ^-0,(r - n) + ^^|T^(r - r.)3/^ 



135 ^5/2 



'' " (B43) 
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This expansion is valid for fixed A 7^ in case |r — r;,| is sufficiently small. Moreover, 
X(r) is an analytic function of (r — rft)^/^ near r = r^. Using this fact, the form of the 
representation (]B43p may be recovered in the usual manner by substitution of an appropriate 
series with undetermined coefficients in differential Eq. (J2]). 
Using the expansion (1B43I) . we have 



^V'^ A, M/2 2 6 + ^2 , , 1^57 + ^2;, ,3/2 

Th' ^n 9 2i/2rf/2 

3105 + 954^2 -llA^ 



F = r'X{r)XAr) - 1 = (-) A{r - nf'' + ^^^f^(r - n) + ^ ^^' + , \ r - r,) 

Vrh/ 3rh ' ' 9 2i/2r?'/2 

-(r-rfe)^ + 0((r-n)5/2). 
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Appendix C: Time Reversal 



(B44) 



This paper has mainly treated geodesic motion in the gravitational ffeld given by met- 
ric ([6]), which was obtained from metric ((51) by time reversal. The purpose of this appendix 
is to study the geodesies of metric (|5]) with X{r) given by Eq. ([2]). As Eq. (|2]) is invariant 
under X i— )■ —X, we can, with no loss in generality, replace X by —X in metric ([5]); that is, 

- ds^ = e-'^i-X^dt" + ^dr^) + —i-Xdt + d(t)f + e'dz^. (CI) 

X j.i J. 

We are interested in the geodesies of Eq. (ICip . If in metric ([6]) we change exp(t) to exp(— t) 
and vice versa, we get metric (JCip . The same transformation turns out to be valid for 
the main equations of this paper that we need for this appendix. That is, in the analysis 
of metric ( ICll) . we ffnd it convenient to deal exclusively with the A > branch of the 
admissible solutions of Eq. (|2]); then, the same equations essentially follow — such as the 
expressions for the invariants Xi and X2 and the equations of motion for the special timelike 
and null geodesies — except that exp(±t) must be everywhere replaced by exp(+t). Consider, 
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for instance, the special timelike geodesies given by 



^ = ''^' 



t , ^2^1/2 



dt _ 

rjX{rj){e' + Clf'\ (C2) 



The solutions are 



dT 

dr 



K (t) = sinh[--|C^|rj(r-ro) + sinh ^(kq)], 



Here, 



= 0o + X(o)(t-to), (C3) 

z{t) = ZQ + K-{tQ)-K~{t). 

it) = |c.le-*/^ 

and K~(to) = ^cT- We note that as t — )■ oo, k^ — )• 0, r — ;• Tend and 2; — )■ ^end, where 

2 



K 



sinh (Kn 



'end — 1^ I Oiiiii V^O J' 

^^end = 2^0 + -ft" (^o)- 



(C5) 



On the other 



It follows that the special timelike geodesies are incomplete in this case [2 
hand, for metric (ICip . 

Xi(o)=4e2Vj, X2(o) = 0, (C6) 

so that as t — )■ 00, Ii — ;• 00. More generally, Xi and X2 are proportional to exp(2t) and 
exp(3t), respectively, in the case of metric (JCip . This is consistent with the incompleteness 
of the special timelike geodesies for t — ?■ 00, since the whole spacetime region becomes 
singular in this limit. Let us note, however, that the special null geodesies are complete, 
since 

t-to = o|C'.|(C-Co), 

0-0o = oX(o)|C',|(C-Co), (C7) 

^-^o = -^(e-*»-e-0, 

so that z{t = 00) = zq + Czexp{—to)/{rj\Cz\)- 

We conclude that the special (timelike and null) geodesies in this case simply wrap in- 
finitely around the axis of symmetry, while z remains forever finite. 
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